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LINEARIZATIONS FOR ROSENBROCK SYSTEM POLYNOMIALS 
AND RATIONAL MATRIX FUNCTIONS 


RAFIKUL ALAM * AND NAMITA BEHERA t 

Abstract. Our aim in this paper is two-fold: First, for computing zeros of a linear time- 
invariant (LTI) system E in state-space form, we introduce a “trimmed structured linearization^’, 
which we refer to as Rosenbrock linearization, of the Rosenbrock system polynomial tS(A) associated 
with S. We also introduce Fiedler-like matrices for <S(A) and describe constructions of Fiedler-like 
pencils for <S(A). We show that the Fiedler-like pencils of »S(A) are Rosenbrock linearizations of the 
system polynomial »S(A). Second, with a view to developing a direct method for solving rational 
eigenproblems, we introduce “linearization” of a rational matrix function. We describe a state-space 
framework for converting a rational matrix function G(A) to an “equivalent” matrix pencil L(A) of 
smallest dimension such that G(A) and L(A) have the same “eigenstructure” and we refer to such a 
pencil L(A) as a “linearization” of G(A). Indeed, by treating G(A) as the transfer function of an LTI 
system S<y in state-space form via state-space realization, we show that the Fiedler-like pencils of 
the Rosenbrock system polynomial associated with Tiq are “linearizations” of G(A) when the system 
Eg is both controllable and observable. 

Key words. LTI system, zeros, poles, rational matrix function, eigenvalue, spectrum, eigenvec¬ 
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1. Introduction. The dynamical behavior of a linear system is often governed 
by a set of linear differential and algebraic equations. We consider a linear time- 
invariant (LTI) system E in state-space form (SSF) given by [M] 

Ex{t) = Ax{t)-i-Bu{t) 
y{t)=Cx{t)+P{i)u{t), 

where A and E are r-by-r matrices with E being nonsingular, P{X) is an n-by-n 
matrix polynomial and B, C are matrices of appropriate dimensions. The system E 
in SSF is a prototype of general LTI systems whose dynamical behaviors are governed 
by matrix polynomials. Indeed, it is shown by Rosenbrock m that a general LTI 
system whose dynamical behavior is governed by a set of equations given by 

T{St)x{t) = U{5t)u{t) and y{t) = V{5t)x(t) + W{5t)u{t) 

can be reduced to an LTI system in SSF such as E by means of strict system equiv¬ 
alence, where r(A), C/(A), U(A), IU(A) are matrix polynomials of appropriate dimen¬ 
sions with T(A) being regular and 6t '■= ^ is time-derivative operator. We, therefore, 
consider LTI systems only in state-space forms. 

The transfer function [14j of the LTI system E is a rational matrix function G(A) 
given by G(A) := P{X) -\- C{XE — A)~^B. The zeros of the transfer function G(A) are 
called the transmission zeros of the LTI system E and are defined as follows mm- 
If diag((?!)i(A)/'!/;i(A), • • • , (^fe(A)/^/>fc(A), 0, • • ■ ,0) is the Smith-McMillan form of G(A) 
then the zeros of the polynomials (t>i{X ),..., (j)k{X) are called the zeros of the rational 
matrix function G(A) [rfl fTOl 1^. 
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The Rosenbrock system polynomial (also referred to as the Rosenbrock system 
matrix) associated with the LTI system E is an (n + r) x (n + r) matrix polynomial 
S{\) given by [TTl [TO] 


S{\) 


[ PW 

C 

B 

{A - XE) 


( 1 . 2 ) 


Note that the LTI system E can be rewritten as S{5t) 


u{t) 

x{t) 



. For an 


ansatz 


u 

X 

we have 5(A) 

u 

X 

gAt ^ 

’ yi.t) ' 
0 

. Hence the ansatz 

u 

X 


produces zero output in the system E when A is an eigenvalue of iS(A) and ^ 

is a corresponding eigenvector. The eigenvalues of iS(A) are called invariant zeros of 
the LTI system E UM- Various zeros such as transmission zeros, invariant zeros 
and decouvlinq zeros of LTI systems play an important role in the analysis of Linear 
Systems 

Computation of zeros of an LTI system is a major task in Linear Systems Theory. 
For the special case when P(A) in (11.11) is a constant matrix polynomial, that is, when 
P{X) := D, the transmission zeros of E can be computed, for example, by matlab 
software command tzero. However, for a general LTI system in SSF such as the LTI 
system E, computation of zeros of the LTI system is a challenging task. We, therefore, 
consider the following problem. 


Problem-I. Compute zeros (transmission zeros, invariant zeros, and decoupling 
zeros) of an LTI system E in SSF. 


We develop a direct method for computing zeros especially transmission zeros 
and invariant zeros of an LTI system in SSF. An obvious direct method for comput¬ 
ing invariant zeros of the LTI system E is to consider a linearization [3111] of the 
Rosenbrock system polynomial 4S(A) and then solve the resulting generalized eigen¬ 
value problem. However, such a linearization based method has serious drawbacks. 
Firstly, if the degree of the polynomial P(A) is m then a linearization of S{\) would 
result in an m{n -|- r) x m(n r) matrix pencil L(A) which would be too large for a 
modest value of m when n and r are large. Secondly, such a linearization L(A) would 
have (spurious) eigenvalues at oo when m > I even when the leading coefficient of 
P(A) is nonsingular. It is well known that infinite eigenvalues of a pencil often pose 
numerical difficulty in computing eigenvalues of the pencil accurately. Thus compu¬ 
tation of eigenvalues of the pencil L(A) with multiple eigenvalues at infinity would 
be inadvisable. Thirdly, although the Rosenbrock system polynomial iS(A) is a struc¬ 
tured matrix polynomial, a linearization L(A) of SiX) may not preserve the structure, 
that is, the pencil L(A) may not be the Rosenbrock system matrix of an LTI system 
in SSF. 

To overcome the drawbacks just outlined, we introduce a trimmed structured lin¬ 
earization of the Rosenbrock system polynomial S{\) which we refer to as Rosenbrock 
linearization of 5(A). More precisely, we introduce Rosenbrock linearization of 5(A) 
which is an {mn + r) x {mn + r) matrix pencil T(A) such that T(A) is the Rosenbrock 
system matrix of an appropriate LTI system in SSF and that T(A) is equivalent to 
5(A) in the sense of system equivalence. We also introduce Fiedler-like pencils for the 
system polynomial 5(A) and show that the Fiedler-like pencils are Rosenbrock lin¬ 
earizations of 5(A). We show that a Rosenbrock linearization of 5(A) deflates infinite 
eigenvalues of 5(A) when the leading coefficient of P(A) is nonsingular. 
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Analysis of transmission zeros of LTI systems naturally leads to zeros of rational 
matrix functions [Humiin]. On the other hand, rational eigenvalue problems arise 
in many applications such as in acoustic emissions of high speed trains, calculations of 
quantum dots, free vibration of plates with elastically attached masses, and vibrations 
of fluid-solid structures, see [Il[Tili[2l[Tll|T3l[3l[2a[2T] and references therein. 
For example, the rational eigenvalue problem (REP) 

k ^ 

G{X)x := —Ax + XBx + X^ ^ = 0, 
arises when a generalized linear eigenproblem is condensed exactly |12) , and the REP 

k 

G{X)x := {K - XM + ^ —^Cj)x = 0, 

j=l ^3 

arises in the study of the vibrations of fluid solid structures, where K = and 
M = M'^ are positive definite and Gj = Gj,j = 1 : fc, have low ranks [13 El [22]. 

Thus computation of eigenvalues and eigenvectors of a rational matrix function 
is an important task which requires development of efficient numerical methods. The 
development of numerical method for solving a rational eigenvalue problem is an 
emerging area of research [1313 [13 E3I13 mils]- We therefore consider the following 
problem. 

Problem-II. Develop a direct method for solving an REP G{X)u = 0, where 
G(A) is a rational matrix function. 

An obvious direct method to solve an REP is to transform the REP to a poly¬ 
nomial eigenvalue problem (PEP) by clearing out the denominators in the rational 
matrix function followed by linearization of the resulting PEP to obtain a generalized 
eigenvalue problem (GEP). Schematically, 

REP —^ PEP —^ GEP —^ SOLUTION. 


A downside of this brute-force “polynomialization” of an REP is that the transforma¬ 
tion from REP to PEP may introduce spurious eigenvalues, which may be difficult 
to detect and remove. Moreover, the transformation from REP to PEP followed by 
linearization may result in a GEP of very large dimension especially when the ratio¬ 
nal matrix function has a large number of poles m- For an illustration of spurious 
eigenvalues, consider the rational matrix function 


G(A):= 


1 

0 


1 • 

X-2 

1 


(1.3) 


Note that G(A) does not have an eigenvalue (oo is not an eigenvalue of G(A)) whereas 
its polynomialization P(A) := (A — 2)G(A) has a multiple eigenvalue at A := 2. 

On the other hand, nonlinear eigensolvers such as Newton method and nonlinear 
Rayleigh-Ritz methods (e.g., nonlinear Arnoldi, rational Krylov, Jacobi-Davidson) 
may be suitable when a few eigenpairs are desired but their convergence analysis is a 
challenging task, see [13 [23 El [23 [24] [HI and references therein. 

Thus a better alternative would be to “linearize” an REP to obtain a GEP of 
least dimension in such a way that avoids a polynomialization of the REP. Such a 
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method based on minimal realization of a rational matrix function is proposed in 
which avoids polynomialization of REP. Indeed, by considering a minimal realization 
of G'(A) of the form 


G(A) = y: X^Aj + C{\E - A)-^B, (1.4) 

3=0 


where A, E, C, B are constant matrices of appropriate dimensions with E being non¬ 
singular, it is shown in HZ] that the eigenvalues and eigenvectors of G'(A) can be 
computed by solving the generalized eigenvalue problem for the pencil 



-^m 



■^m—1 ■^m — 2 ' ’ ' -^0 
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-E _ 


B 

A 


(1.5) 


where the void entries represent zero entries. The pencil C(A) is referred to as a 
companion linearization of G'(A) in [17] and can be thought of as an extension of the 
companion linearization 


Gp(X) := X 

Am 0 

0 In 

... O' 

+ 

Am — l Am— 2 
-In 0 

Aq 

0 


_ 0 

1 

o 


0 

-In 0 


of the matrix polynomial P(A) := the case of G(A). The pencil Cp{X) 

is also known as the first companion form of P(A) [7]. The size of the pencil C(A) 
is usually much smaller than that of a pencil obtained by polynomialization of G(A) 
followed by linearization specially when the coefficient matrices of G(A) have low 
ranks [TT] . 

The pencil C(A) has been referred to as “companion linearization” of G(A) in [17]. 
So it is natural to ask: Is the pencil C(A) a linearization of G{X) ? If yes, then in what 
sense does C{X) define a linearization of G{X)7 It is well known that the companion 
form Gp(A) is a linearization of P(A) in the sense of unimodular equivalence [7l ITT]. 

With a view to addressing Problem-II as well as the questions raised above, we 
consider yet another problem. 

Problem-III. Let diag((()i(A)/'!/'i(A), • • • , (/)fc(A)/'0fe(A), 0, • • • ,0) be the Smith- 
McMillan form of G(A). Does there exist a matrix pencil L(A) of smallest dimension 
such that diag(/p, (j)i{X), - ■ ■ , (j)k{X),0, ■ ■ ■ ,0) is the Smith form of L(A) for some ap¬ 
propriate p7 Describe a framework for construction of such a pencil when it exists. 

The pencil L(A), when exists, can be considered as a “linearization” of G(A) 
as it preserves the spectrum (zeros) of G(A) and the partial multiplicities of the 
eigenvalues. Observe that L(A) exists and is a linearization of G(A) when G(A) is a 
matrix polynomial. Consequently, the linearization in Problem-III, when exists, is 
an extension of linearization of a matrix polynomial to the case of a rational matrix 
function. 
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We provide a complete solution of Problem-III. We show that computation 
of eigenvalues and eigenvectors of a rational matrix function is closely related to 
computation of transmission zeros of an LTI system. We therefore embed an REP 
into an appropriate LTI system in SSF. More precisely, we treat a rational matrix 
function G'(A) as the transfer function of an LTI system E in SSF via state-space 
realization [14]. In such a case, the zeros of G'(A) become the transmission zeros of 
the LTI system S. Thus, we compute transmission zeros of the LTI system E from 
which we recover the eigenvalues of G{X). We show that a Rosenbrock linearization 
L(A) of the LTI system E is a linearization of G'(A) in the sense described in Problem- 
III when the LTI system E is controllable as well as observable. In particular, we show 
that the companion form C(A) in (11.51) is indeed a linearization of G(A) in the sense 
described in Problem-III when the realization of G(A) in (II.4p is minimal. Thus, 
schematically, our strategy for solution of Problem-II can be summed as follows: 

Rational matrix function —> State-space realization —> Linearization —> Solution. 

Finally, we mention that the set of transmission zeros of an LTI system E forms 
a subset of the set of invariant zeros of the LTI system and the transmission zeros 
coincide with the invariant zeros when the LTI system E is controllable as well as 
observable. Thus, Problem-II via state-space realization is subsumed by Problem- 
1. There are efficient methods for computing a (minimal) state-space realization of a 
rational matrix function [uiinKii]- For example, the MATLAB command minreal can 
be used to compute a minimal realization. We therefore focus on solving Problem-I 
and address the issues in Problem-III by assuming a minimal state-space realization. 

The rest of the paper is organized as follows. Section 2 describes spectra and poles 
of rational matrix functions. Section 3 introduces Fiedler matrices associated with a 
Rosenbrock system matrix and describes construction of Fiedler pencils. Section 4 
introduces Rosenbrock linearizations for Rosenbrock system matrix and shows that 
Fiedler pencils are indeed Rosenbrock linearizations for Rosenbrock system polyno¬ 
mials. Finally, Section 5 introduces linearization of a rational matrix function, which 
extends the notion of linearization of a matrix polynomial to the case of a rational 
matrix function, and describes construction of such a linearization via minimal state- 
space realization of the rational matrix function. 

Notation. We denote by C[A] and C(A), respectively, the polynomial ring and the 
field of rational functions of the form p{X)/q{X), where p(A) and q{X) are scalar poly¬ 
nomials in C[A]. Further, we denote by C[A]™^" and C(A)'"^"', respectively, the 

vector spaces (over C) of m-hy-n matrices, matrix polynomials and rational matrix 
functions. We denote the Smith form [H] of ^(A) G C[A]'”^" by SF(X(A)) and the 
Smith-McMillan form [T3j of y(A) G C(A)'"^" by SM(y(A)). The normal rank [2] of 
X(A) G C(A)'”^" is denoted by nrk(X) and is given by nrk(X) := max^ rank(X(A)) 
where the maximum is taken over all A G C which are not poles of the entries of ^(A). 
We denote the j-th column of the n x n identity matrix /„ by Cj and the transpose of 
a matrix A by . We denote the Kronecker product of matrices A and B hy A^ B. 

2. Eigenvalues, zeros and poles. Let P{X) G C[A]™^”. If nrk(X) = min(m, n) 
then P{X) is said to be regular, otherwise P(X) is said to be singular. A matrix 
polynomial U{X) G C[A]"'^" is said to be unimodular if det(C/(A)), the determinant 
of U{X), is a nonzero constant for all A G C. Suppose that the Smith form [14] of 
P(A) is given by SF(P(A)) = diag(^i(A),... ,(()fe(A),0,... ,0), where (fiiX), ... ,(()fc(A) 
are invariant polynomials (unique monic polynomials such that (fi divides 4>i+i for 
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i = 1,..., fc — 1) of P(A) and k = nrk(P). We define the zero polynomial (j)p{\) of 
P(A) and the spectrum (also called the zeros) Sp(P) of P(A) by 


k 

</,p(A) := ]^<))j(A) and Sp(P) := {A e C : (?ip(A) = 0}. (2.1) 

i=i 


A complex number A G C is said to be an eigenvalue of P(A) if rank(P(A)) < nrk(P). 
It follows that /i is an eigenvalue of P(A) if and only if (()p(/i) = 0. Thus the spectrum 
Sp(P) consists of (finite) eigenvalues of P(A). We mention that the spectrum Sp(P) 
as defined above does not include infinite eigenvalues of P(A). 

Now consider the LTI system S given in (HID and the associated Rosenbrock 
system polynomial 5(A) given in (II.2|) . Then the spectrum Sp(5) of 5(A) is called the 
invariant zeros of the LTI system E. If A G Sp(5) then a corresponding eigenvector 
with u G C" and cc G C*" is called an invariant zero direction and, the 
components u and x are, respectively, called input direction and state direction of the 
LTI system E. Further, the spectrum of the matrix pencils 


M(A) := [A - \E B] and N{X) := 


A - XE 
C 


( 2 . 2 ) 


are, respectively, called the input decoupling zeros and output decoupling zeros of the 
LTI system E. We refer to mumii] for more on zeros of an LTI system. 

Next, let G'(A) G C(A)’”^" be a rational matrix function. Suppose that the 
Smith-McMillan form SM(G(A)) of G(A) is given by [TTl [T^ 

SM(G(A)) = diag (^1^, • • • , . (2.3) 

where the scalar polynomials ^i(A) and '0i(A) are monic, are pairwise coprime and, 
4>i{,X) divides (jii+i (A) and '0i+i(A) divides '0i(A), for z = 1,2,..., /c — 1. The polyno¬ 
mials (/'i(A),..., 4>k{X) and '0i(A),. ■ •, V'fc(A) are called invariant zero polynomials and 
invariant pole polynomials of G(A), respectively. Define the zero polynomial 4>gW 
and the pole polynomial tpoW of G(A) by 

k k 

<(>g(A) := 0j(A) and V'g(A) := nV'j(A). (2.4) 

i=i 1=1 

Definition 2.1 (Zeros and poles [HI). Let G(A) G C(A)™^”. Let (pcW and 
ipcW be the zero and pole polynomials o/G(A), respectively. A complex number X is 
said to be a zero of G(A) if 4 >gW = 0. A complex number X is said to be a pole of 
G(A) iffjG{X) = 0. 

We mention that G(A) is said to have a zero (pole) at A = oo if G(l/A) has a zero 
(pole) at A = 0. We, however, consider only finite zeros and poles of G(A). It follows 
that a complex number /r is a pole of G(A) if and only if ||G(A )||2 —>■ oo as A ^ We 
denote the set of poles of G(A) by Poles(G) and define the spectrum Sp(G) of G(A) 
to be the set of zeros of G(A). In other words, we set 

Poles(G) := {A G C : V'g(A) = 0} and Sp(G) := {A G C : <^g(A) = 0}. (2.5) 
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The spectrum of the transfer function of an LTI system is called the transmission zeros 
of the LTI system M- Thus, the transmission zeros of the LTI system E in CH) are 
given by the spectrum of the transfer function G(A) := P(A) + C{XE — A)~^B. 

Definition 2.2 (Eigenvalue and eigenpole). Let G(A) € C(A)'"^". A complex 
number Xq € C is said to be an eigenvalue of G{X) i/rank(G(Ao)) < nrk(G). We 
denote t/ie eigenspectrum (set of eigenvalues) of G{X) &?/Eig(G). 

A complex number Aq is said to he an eigenpole of G(A) if Aq is a pole of G(A) 
and there exists v{X) G C"[A] such that ?;(Ao) ^ 0 and limA->.A(, G(A)?;(A) = 0. We 
denote the set of eigenpoles of G{X) by Eip(G). 

We mention that G(A) is said to have an eigenvalue (eigenpole) at A = oo if 
G(l/A) has an eigenvalue (eigenpole) at A = 0. Eor example, G(A) in (11.31) has no 
eigenvalue but has an eigenpole at A = 2. On the other hand, G(A) := diag(A, 1/A) 

has no eigenvalue but has eigenpoles at 0 and oo. We, however, consider only finite 

eigenvalues and eigenpoles of rational matrix functions. 

Theorem 2.3. Let G(A) e C(A)’”''”. Then Sp(G) = Eig(G) U Eip(G) and 

Eig(G) = {A e C : (/g(A) = 0 and i/g(A) + 0}, 

Eip(G) = {A G C : /'g(A) = 0 and V'g(A) = 0}, 

where 4 >g{X) and '0g(A) are zero and pole polynomials o/G(A), respectively. 

Proof Suppose that SM(G(A)) = diag ((()i(A)/i/>i(A),..., (()fe(A)/'0fc(A), 0,..., 0). 
Then (/g(A) = 0^=1 ^g(A) = Suppose that Aq G Eig(G). 

Then we have rank(G(Ao)) < k = nrk(G). Consequently, (/i(Ao) = 0 for some i and 
hence </>g(Ao) = 0 . Since cfi and ifi are coprime, it follows that ipi^Xo) 0 for i = 1 ■. k 
and hence V'g(Ao) ^ 0 . Conversely, if 4>g{Xo) = 0 and 'i/'g(Ao) ^ 0 then />i(Ao) for 
some i and hence rank(G(Ao)) < k = nrk(G) showing that Aq G Eig(G). 

Next, let Aq G Eip(G). Then there exists v{X) G C”[A] such that u(Ao) 7 ^ 0 and 
limA^Ao G(A)u(A) = 0. By Smith-McMillan form, G(A) = C/(A) SM(G(A)) 1/(A) for 
some unimodular matrix polynomials 17(A) and 17(A). Hence 

U-'^(X)G{X)v(X) = SM(G(A))17(A)u(A). (2.6) 


Setting a;(A) := 17(A)u(A) it follows that a;(A) — )• a;(Ao) = 17(Ao)u(Ao) ^ 0. Hence by 
(1^ we have ^‘|^j xi(A) —>• 0 for i = 1 : fc. Since a;(Ao) 7^ 0, we have a:i(Ao) ^ 0 for 
some i. This shows that as A —7 Aq, we have 




1 

Xi{X) '(/'j(A) 


Xi{X) —7 0. 


Since (pi and ipt are coprime, 'ipi(Xo) 0 and (j)i{X) (/>z(Ao) = 0 for some i. Hence 

</>g(Ao) = 0. Also V’g(Ao) = 0 since Aq is a pole. 

Conversely, suppose that 4>g{Xo) = 0 and V'g(Ao) = 0. Define u{X) := 17“^(A)ei, 
where 17(A) and V{X) are as in (j2.6l) . Then 


G(A)u(A) = t/(A)SM(G(A))ei = 17(A)44^ei ^ 0 as A ^ Aq, 

V'i(A) 


which shows that Aq G Eip(G). Hence the proof. □ 

Thus the spectrum of a rational matrix function G(A) consists of eigenvalues and 
eigenpoles. We define partial multiplicities of zeros (eigenvalues and eigenpoles) of a 
rational matrix function as follows. 

Definition 2.4 (Partial multiplicities of zeros). Let G(A) G C(A)"*^" and 
(pi{X),..., (j)k{X) be invariant zero polynomials ofG{X). Let Xq G Sp(G). Then (Pg{Xo) = 
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0 and = {X — Xo)'^*di(X) with di{XQ) ^ 0 and'ji > 0 fori = 1 : k. The index tuple 
Ind0(Ao,G') := (71 ,... ,7fc) is called the multiplicity index of G at Ag and satisfies the 
condition 0 < 71 < 72 < • • • < 7fe. The nonzero components in Ind^(Ao, G) are called 
the partial multiplicities of Ag as a zero of G{X). The factors (A — Ag)'*'* with 7i 0 
are called the elementary divisors of G{X) at Ag. The algebraic multiplicity of Xq is 
defined by 

0'<j>{Xo) := 7i + 72 + • • • + 7 fe = multiplicity of Ag as a root of (j)c[X). 

If a^{Xo) = 1 then Ag is called a simple zero of G{X). 


Similarly, we define partial multiplicities of poles of a rational matrix function as 
follows. 

Definition 2.5 (Partial multiplicities of poles). Let G{X) G C(A)™^" and 
'if>i{X),... ,ipk{X) be invariant pole polynomials of G{X). Let Ag G Poles(G'). Then 
= 0 and '0i(A) = (A — Ag)“’gi(A) with (?i(Ag) 0 and ai > 0 for i = 1 , 2 ,... ,k. 

The index tuple Ind^(Ag, G) := {ak, a^-i . ■., oli) is called the multiplicity index of G 
at Ag and satisfies the condition ak < ctk-i < ■ • • < cii. The nonzero components in 
Ind. 0 (Ag,G) are called the partial multiplicities of Xq as a pole of G{X). The factors 
(A — Ag)“* with ai ^ 0 are called the elementary divisors of G{X) at the pole Ag. The 
algebraic multiplicity of Ag is defined by 

a-ipiXo) := Q!i + a2 + ■ • • + cife = multiplicity of Xq as a root ofipciX). 

If cLijiiXo) = 1 then Ag is called a simple pole of G{X). 


We remark that if Ag G Eip(G) then Ind,^(Ag,G) yf 0 and Ind^(Ag,G) 7 ^ 0. In 
such a case, Ag has two multiplicity indices, namely, the multiplicity index as a zero 
and the multiplicity index as a pole. In general Ind0(Ag,G) yf Ind.0(Ag,G). 


Example 2.6. Consider G(A) := 


A(A-2)^ 


A-2 

A 


. Then X = 2 is an eigenpole 


with Ind0(2, G) = (0,1) and Ind^(2, G) = (0, 2). Hence Ind0(2, G) ^ Ind^(2, G). 


3. Fiedler pencils for Rosenbrock system polynomial. For the rest of the 
paper, we consider the linear time-invariant (LTI) system E given in (11.11) . that is. 


Ex{t) = Ax{t) + Bu{t) 
y{t) = Gx{t) + P{-^)u{t), 


where E G is nonsingular and P{X) G C"^"[A]. Also, unless stated otherwise. 


5(A) 


[ PW 

C 

B 

1 


and G(A) := P(A) -k G{XE - A)-^B 


(3.1) 


will denote the Rosenbrock system polynomial and the transfer function of the LTI 
system E, respectively. We also refer to 5(A) as the Rosenbrock system matrix of the 
LTI system E. We define the degree of the Rosenbrock system polynomial 5(A) to 
be the degree of the matrix polynomial P{X). For the rest of the paper, we assume 
that the degree of P{X) is m and P{X) is given by 


P(A) :=^A^A, 
0=0 
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with Am yf 0. 


(3.2) 









We now introduce Fiedler pencils of 5(A) which are potential candidates for 
“trimmed structured linearizations” of 5(A) and describe their constructions. We 
mention that although the pencil A — XE in 5(A) is assumed to be regular and E to 
be nonsingular, which are important for the associated LTI state-space system, these 
assumptions are irrelevant for construction of Fiedler pencils of 5(A). The construc¬ 
tion of a Fiedler pencil described below holds good even when 5(A) and P{X) are 
singular polynomials and A — XE is a singular pencil. 

For a ready reference, we briefly describe Fiedler matrices associated with P{X) 
and refer to [B] for further details. We closely follow the notational conventions used 
in . The nm x nm matrices Mq ,, Mm given by 


Mn := 


{m —l)n 


4 I Mm ■ — 

— Ar 


d-('m—l)n 


(3.3) 


and 


M, := 




-A, 

In 


In 

0 


hi- 




i = 1,..., TO — 1, (3.4) 


are called Fiedler matrices associated with P{X). The Fiedler matrices Mi and Mj 
commute when \i — j\ > 1, that is, MiMj = MjMi when \i — j\ > 1. Further, for 
i = l,...,TO — 1, each Mi is invertible and 


M-i 


I(m. — i — 1 


0 

In 


In 

Ai 


I{i — l)n 


(3.5) 


Let a : {0,1,..., to — 1} —>• {1,2,..., to} be a bijection. Then M„ denotes the 
product of Fiedler matrices given by 


Ma := Mrj-l(^i)Mcr-i(2) ' ' ' 

Note that a(i) describes the position of the factor Mi in the product Mo-, that is, 
a{i) = j means that Mi is the j-th factor in the product. For convenience, we set 
Mfh := Inm, where 0 is the empty set. The nm x nm pencil Lc{X) given by 


L,{X) := XMm - M, 


is called the Fiedler pencil of P{X) associated with a, see [6l[I]. 

We now introduce Fiedler-like matrices associated with the Rosenbrock system 
polynomial 5(A) and construct Fiedler-like pencils for 5(A). Recall from (13.11) and 
(1321) that 5(A) is an (n -|- r) x (n -|- r) matrix polynomial of degree to. 

Definition 3.1 (Fiedler matrices). Define the {nm -|- r) x {nm + r) matrices 
Mo,...,Mm by 


and 


Mo := 


Mo 

Cm ^ C 

M — 

■ Mm 

0 

-e^(S)B 

-A 

, -l-viirn • 

0 

-E 


(3.6) 


M, 

0 

0 

Ir \ 


9 


i = 1,..., TO — 1, 


(3.7) 






















where Mi,i = 0,1,... ,m are Fiedler matrices associated with P{X)- We refer to the 
matrices M^, i = 0,1,... ,m, as the Fiedler matrices associated with S{X). 

Observe that unlike in the case of Fiedler matrices associated with P{X), the 
Fiedler matrices Mq and do not commute, that is, MoM^ ^ M^Mq. However, 
Mi commutes with Mj when \i — j\ > 1 and \i — j\ ^ rn, that is, 

MiMj = MjMi for \i — j\ > 1 , except for Mm and Mq. (3-8) 


For i = 1 : TO — 1, each M^ is invertible and 



, where ^ 


is given in (USD. 


Hence we have 


(3.9) 


Mj ^M^- ^ = M^- ^Mj ^ for \i — j\ > 1, except for Mq ^ and M^^. (3.10) 

Definition 3.2 (Fiedler pencil). Let a : {0,1,..., to — 1} —>■ {1,2,..., to} be a 
bijection. Then the (mn + r) x (mn + r) matrix pencil Lct(A) given by 


Lct(A) := AMm — M^- 1 (i)M£,- 1 ( 2 ) • • • MCT-l(m) = AMm — Mct (3.11) 

is said to be the Fiedler pencil of the Rosenbrock system polynomial S{X) associated 
with a. 

If <7 and T are two bijections from {0,1,..., to — 1} to {1, 2,..., to} then because 
of commutation relation we may have M^ = M,-. Consequently, a Fiedler pencil may 
be associated with more than one bijection. 

Definition 3.3. Let a, t : (0,1,..., to — 1} {1, 2,..., to} be bijections. Then 

a is said to be equivalent to t (denoted as a ^ t) ifM.fj = M,-. 

Example 3.4. Consider G(A) := A‘^H 4 + • • ■ + Hg + C{XE — A) and the 
Fiedler pencil Lct(A) := AM 4 — M 1 M 0 M 2 M 3 with a~^ = (1,0, 2,3). Then 




—A^ 

In 

0 

0 

0 


■ AI 4 



—A2 

0 

In 

0 

0 

U(A) = A 

hn 


- 

-Ai 

0 

0 

— ^0 

-c 



-E 


In 

0 

0 

0 

0 



0 

0 

0 

-B 

-A 


Considering Lt-(A) := AM 4 — M 2 M 0 M 1 M 3 with t ^ = (2, 0,1,3), we have 


A4 



— ^3 

In 

0 

0 

0 

In 



— A2 

0 

-All 

In 

0 

In 


— 

In 

0 

0 

0 

0 

In 



0 

0 

—aIo 

0 

-c 


-E 


0 

0 

-B 

0 

-A 


For the bijection S given by S ^ = (0,2, 3,1), we have Lt-(A) = L 5 (A). This shows that 
S'^T.m 

Definition 3.5. m Let a : {0,1,..., to — 1} ^ (1, 2,..., to} be a bijection. 

(1) For d = 0,... ,m — 2, we say that a has a consecution at d if a(d) < a(d+l) 
and a has an inversion at d if a{d) > a(d + 1). 
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(2) The tuple CISS{a) := (ci, ii, C 2 , 12 , • ■ • ,ci,ii) is called the consecution-inversion 
structure sequence of a, where a has Ci consecutive consecutions 0 , 1 ,..., Ci — 
1 ; ii consecutive inversions oi ci, ci + 1 ,..., ci + zi — 1 and so on, up to ii 
inversions at m — 1 — ii,... ,m — 2. 

(3) We denote the total number of consecutions and inversions in a by c(a) and 

i i 

i{cr), respectively. Then c(cr) = cj, i(a) = ij, and c{a)-\-i{a) = m—1. 

i=i 

Said differently, if cr : {0,1,..., m — 1} —>■ {1,2,..., m} is a bijection then cr has 
a consecution at d if and only if is to the left of M^+i in M^., while a has an 
inversion at d if and only if is to the right of M^+i in M^. Further, we have 


MiMj = MjMi 


MjMi for i,j € a. 


We now show that a Fiedler pencil of 5(A) can be constructed from a Fiedler 
pencil P{X) and vice-versa. In fact, there is a bijection from the set of Fiedler pencils 
of 5(A) to the set of Fiedler pencils of P(A). 

Theorem 3.6. Let a : (0,1,... , m — 1} ^ {1, 2,..., m} be a bijection. Let L„{X) 
and Lcr(A) be the Fiedler pencils of P{\) andS(\), respectively, associated with a, that 
is, Tcr(A) := XMjn — and Lcr(A) := AM^ — Mo-. If = (erf 0 , erf ^) for some 
bijections cri and a 2 , then 


L.(A) 


r La{X) 

(em ® C) 

[ (e)f 0 B)M,, 

{A - XE) 


Further, if CISS{a) = (ci, ii,..., c;, f) then 


and 


L.(A) 


\ L,{X) 

Cm C 


A-XE 


if Cl > 0 


U(A) 


LAX) 

C(m—ii) ^ ^ 

ef, (8)5 

A-XE 


if Cl = 0. 


Thus the map Fiedler(P) —^ Fiedler(5), XMm — 1 —>■ AM^ — is a bijection, 

where Fiedler(P) and Fiedler(5), respectively, denote the set of Fiedler pencils of P(X) 
and 5(A). 

Proof. We have Lct(A) = AM^ — M^- = AM™ — McriMoMo .2 


' Mm 

0 

_ 

■ Men 

0 

\ Mo 

Cm C 

0 

■ Mm 

-E 

0 

0 

Ma, 

Ir \ 

^MqM 

[ -e 

1 -M,nl 

-A 

[Cm (8> C) ] 

0 

-E 

[ (-e^ (8> 5)M,^2 



LAX) 

Men Am <8 C) 

_ (ef, (g) 5)M,^2 

1 


Merj 

0 

0 



Now suppose that ClSS(er) = (ci, zi,..., c;, f). 

Case / : Suppose that ci > 0. Then by commutativity relation we have = 
Men MqMi ■ • • Mci with Cl -I- 1 G cri. Thus M^r = MoMct 2 , where = Mi • • • M^. 
Hence 


■ Men 


A 

\ Mo 

Cm ® C 

M, 

0 ^ 

1 _ 

. -e 

72 

in®B 

MaA- 

-A 

Cm^C) ' 

. (-e^ <8 5)M,^3 

J 
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Since j G tJi implies that j > ci + 1, we have 


. This shows that 


(em ® In) = In ^nd (— 6^ ® C) = —Cm ® C. Next, we have 


(e™ ® Iu)Mi = (e^ (g) In) 


(e^ (g) /„)MiM 2 = (e^_i (g In) 


^ (m — 2)n 


‘(m—3)n 


In 

In 0 

— A 2 In 

In 0 


= (e^_i (g) J„), 


~ i^m-2 ® ^n) 


and so on. Thus (e^ (g /ji)MiM 2 • • -Md = (eJn-ci ® In)- Hence (e^ g) 
(e^_c g In) and (—g B)M„^ = —{e^_ g B). Consequently, we have 


r -La{X) 

Cm ® C 

_ e(m-c,) ® B 

A-XE 


L,^(A) = AM™ -M,^ = 


Case II : Suppose that c\ = 0. Then a has ii inversions at 0. Hence by commu¬ 
tativity relations we have = Mi^ • ■ -MiMoMo-j =: McriMoMcr 2 with ii -I- 1 £ (T 2 - 
Hence 


\ Mn, 



Ir \ 


\ Mo 

Cm ® C 

_ -e™gH 

-A 


Ma-2 



Ir J 



A7o-i(-e™ g C) 

[ {-ein®B)M„, 

-A 


Since j £ U 2 implies that j > *i -I- 1, we have 
(e^ 

I{m—2)r 





Ii\n 


. This shows 


that (e™ g /„)Mct 2 = e™ g Hence g H) = —g B. Next, we have 


(e™ g In) — 


ll2Mi(em g J„) — 


I{m—3)r 


— Hi In 
In 0 

— H 2 In 
In 0 


(c™ g In) — (e™ —1 g In)j 


(e ™—1 g In) — (e™ —2 g In)- 


Thus Mij • • •M 2 Mi(e™ g /„) = {cm-n ® In)- Hence Mcr^(e™ g In) = (e(™_q) ® In) 
and e™ g C) = ® C)- Consequently, we have 

L^(A) = AM™ -M^ = 


\ -Lo{X) 

Cm—ii ^ C 

_ ein®B 

A-XE 


Recall that for each i,j £ a, we have MiMj = MjMi M^Mj = Mj-M^. Hence it 
follows that #(Fiedler(P)) = #(Fiedler(5)). This completes the proof. □ 

The companion pencil C(A) given in (11.511 is in fact a Fiedler pencil of 5(A). 
Proposition 3.7. Let C(A) be the eompanion pencil given in / H..5I) . Then 


C(A) = AM™ -] 




IlMo — Lcr(A), 
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where a~^ = (m — 1 ,..., 2,1, 0). 

Proof. Consider the Fiedler pencil L£r(A) with a~^ = {m — 1,... ,2,1, 0). Then it 
follows that ClSS(cr) = {0,m — 1). Hence by Theorem 13.61 we have 


XMm — Mm-1 ■ ■ ■ MiMq 

ei®C 


\ Cp{X) 

Cl ® C 


{A - XE) 



-1 

1 


where C'p(A) is the companion pencil of P(A) given in (|1.6I1 . The last equality follows 
from the fact that C'p(A) = XM^ — Mm-i • ■ • MiMq, see [6]. This shows that Lct(A) = 
C(A) as desired. □ 

The pencil C(A) is derived in |17) and is referred to as companion linearization 
of G'(A). We refer to C(A) as the hrst companion form or the first companion pencil 
of 5(A). We also refer to C(A) as the first companion pencil of the transfer function 
G'(A) associated with 5(A). 

We define the second companion form of 5(A), which we denote by C2(A), by 


C2(A) — AMm — MqMi ■ • ■ — L(t(A), 


where = (0,l,...,m — 1). Since ClSS(cr) = (m — 1,0) and C'2(A) = XMm — 
Mq ■ ■ ■ is the second companion form of the polynomial P{X), see [6], by 

Theorem 13.61 we have 


r c 2 (a) 

Cm 0 C 

e ( (g) B 

-< 

1 





^m—1 In 


In 



— Am-2 0 





i ■■ • In 


In 



O 

O 

o 

1 

-C 


-E \ 


-B 

-A 


(3.12) 


We refer to C2(A) as the second companion form of 5(A) or the transfer function C?(A) 
associated with 5(A). 

It is well known that a matrix polynomial admits block pentadiagonal Fiedler 
pencils [6]. The next result characterizes block pentadiagonal Fiedler pencils of Rosen- 
brock system polynomials. 

Theorem 3.8. Let a : {0,1,..., m — 1} {1,2 ,...,to} be a bijection. Let 

Lct(A) := XMm — Mu and Lcr(A) := AM^ — Mg-, respectively, be the Fiedler pencil 
of P{X) and S{X) associated with a. Suppose that ClSSljr) = {ci,ii,... ,ci,ii). Then 
Lcr(A) is block pentadiagonal if and only if Lu{X) is block pentadiagonal and ci < 1 
and ii < 1. In .such a case, we have 


U(A) 


[ LaW 

Cm ® C 

. e^_i(S)B 

1 


*/ Cl > 0 


and 


L.(A) = 


r LaW 

Cm—1 ^ C 

_ e^(S>B 

1 


if Cl = 0. 
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Proof. The proof is immediate by Theorem l3.6l Indeed, for pentadiagonal Lcr(A), 
we must have Ci < 1 and ii < 1. If Ci > 0 then Ci = I and ® B = 0 B. 

So by the Theorem 13.61 we have 


r LaW 

Cm C 


r La{X) 

Cm ^ C 

. ® B 

-< 

1 



{A - XE) 


is block pentadiagonal if and only if To-(A) is block pentadiagonal. Similar argument 
holds when ci = 0 and ii = I. □ 

Let O = MiMz • • • be the product of odd Mi factors and £ = M 2 M 4 • • • be the 
product of the even Mi factors, excluding Mq and Mm- Then it is shown in [6] that 
the product of O, £ and Mg in any order gives a pentadiagonal Fiedler pencil of T’(A). 

Corollary 3.9. Let he the product of odd factors and be the product 
of even factors, excluding Mg and M^. Then Lo.(A) = AM^ — Mo-iMgMcra and 
L5(A) = AMm — MoM^sMcri are block pentadiagonal pencils of S{X). However, the 
pencils L(A) = AM^ — M^aMo-iMg and L(A) = AM^ — MoMcnM^j are not block 
pentadiagonal. 

Proof Consider the pencil Lcr(A) = AM™ — M^jMgMcrs. Then ClSS(cr) = 
(ci, ii,...) = (0,1,...). Hence by Theorem 13.61 we have 


\ 

Cm—ii ^ C 


r LaW 

Cm —1 ^ C 

_ e^(£)B 

1 

-< 

1 



-1 

-< 

1 


where To-(A) = \Mm — It is well known that the product of 

Mq in any order gives a block pentadiagonal matrix (i, Example 3.2). Hence Lcr(A) = 
XMm — Mcr is block pentadiagonal and consequently Lcr(A) is block pentadiagonal. 
For the pencil L5;(A) we have CISS(ct) = ...) = (1,1,...). Hence by Theorem 

13.61 L3:(>) is block pentadiagonal. For L(A) = AM^ — Mg, we have ClSS(cr) = 

(cijii,...) = (2,...). By Theorem 13.61 L(A) has more than two subdiagonal blocks 
and hence is not block pentadiagonal. Similar argument holds for the pencil L(A) = 
AM„ - MgMo.jMcr2. □ 

Example 3.10. Let G(A) := AqX^ + • • • + XAi + Hg + C{XE — A)~^B. Also, let 

= (1,3,5) and = (2,4). We have 



Case I: Consider the pencil Lct(A) := AMg — McriMgMCTa- Then we have 


MgM,,, 


In 


^4 


In 0 


— A 2 In 


In 0 


— Aq 

-C 

-B 

-A 
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and hence 



-A14 

In 

0 

0 

0 

0 

In 

0 

0 

0 

0 

0 

0 

0 

-A13 

0 

— ^2 

In 

0 

0 

0 

In 

0 

0 

0 

0 

0 

0 

0 

0 

-Ai 

0 

—^0 

-c 

0 

0 

0 

In 

0 

0 

0 

0 

0 

0 

0 

0 

-B 

-A 


from which it follows that Lcr(A) is block pentadiagonal. Note that the number of 
consecutions at 0 in is zero and the number of inversion at 0 in is 1. 

Case II: Consider the pencil Lo-(A) := AMg — MoMo-jMo-i. Then 



—^5 

In 

0 

0 

0 

0 

0 


-A4 

0 

—A3 

In 

0 

0 

0 


In 

0 

0 

0 

0 

0 

0 


0 

0 

—A 2 

0 

-Ai 

In 

0 


0 

0 

In 

0 

0 

0 

0 


0 

0 

0 

0 

— Aq 

0 

-C 


0 

0 

0 

0 

-B 

0 

-A 

is block pentadiagonal. Thus Lcr(A) 

is a 

block pentadiagonal pencil. 

number of consecutions at 0 

in Mct is 

1, i.e.. Cl = 

= 1. 




Case III: Consider the pencil Lct(A) 

= AM 

g — MoMo-i 

Mo, 

. The: 

pentadiagonal, since the number of consecutions at 0 

in Mo 

is 2. Inde 


—A 5 

— A 4 

In 

0 

0 

0 

0 


In 

0 

0 

0 

0 

0 

0 


0 

—^3 

0 

— A 2 

In 

0 

0 


0 

In 

0 

0 

0 

0 

0 


0 

0 

0 

-Ai 

0 

In 

0 


0 

0 

0 

—Ao 

0 

0 

-c 


0 

0 

0 

-B 

0 

0 

-A 


which is not block pentadiagonal. Note however that ■= which is 

the (1,1) block matrix in Mo-, is pentadiagonal. Similarly, Mo- = Mo-jMctiMo is also 
not block pentadiagonal, since the number of inversions at 0 in Mcr is 2. ■ 

As in the case of Fiedler pencils of matrix polynomials, Fiedler pencils of Rosen- 
brock system polynomials can be constructed algorithmically. 

Theorem 3.11. Let a : {0,1,..., m — 1} — >■ {l,2,...,m} be a bijection. Let 
Lcr(A) = AMm — Mcr he the Fiedler pencil of system matrix 5(A) associated with a. 
Consider the matrices Wq, Wi, ..., Wm -2 constructed by Algorithm[J\ Then we have 
Mcr = Wm-2- 
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Algorithm 1 Construction of for Lo-(A) := 


m I 1 

T , c,^^ EamJ c \ . 

Input: 5(A) = 


[ B \ A- \E 

{ 1 , 2 ,. 

Output: Mcr 

if cr has a consecution at 0 then 

— Ai In 0 

Wo := -Ao 0 -C 

-B 0 -A 


and a bijection a : {0,1 ,... ,m — 1} 



■ -All 

— Aq 

-C ' 

:= 

In 

0 

0 


0 

-B 

-A 


end if 

for z = 1 : m — 2 do 

if CT has a consecution at z then 

_ Ai-\-l In 0 0 

[ Wi_i(:,l) 0 W,_i(:,2 :t + l) W,_i(:,z + 2) 

else 


W, := 

end if 
end for 

Mct := Wm-2 


-A,+i Wi_i(l,:) 

In 0 

0 W,_i(2:i + 1,:) 

0 W,_i(z + 2,:) 


Proof. We prove the result by induction on the degree m. For m = 2 the proof 
is obvious, since for m = 2 there are only three Fiedler factors and two possibility 
for Mcr, either M^- = MqMi if a has a consecution at 0 or M^- = MiMq if cr has an 
inversion at 0. Now 



■ -All 

In 

0 

MIqMIi = 

— ^0 

0 

-c 


-B 

0 

-A 


' -All 

—aIq 

-c ■ 

In 

0 

0 

0 

-B 

-A 


This proves the case for to = 2. Suppose that the result is true for to — 1 > 2. We 
show that it is true for to and the bijection cr : (0,1,..., to — 1} —>■ {1, 2,..., to}. Note 


satisfies 


li = diag(/„. Mi), for z = 0 : TO — 2, 



m 

I] A^ Aj 

C 

— 

i=o 



B 

A-XE _ 

0 

TO — 2, 



where Mi are the {n{m — 1) + r) x {n{m — 1) + r) Fiedler matrices associated with 

c 1 


S(A) = to 


B \ A-\E 
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Case I. If a has a consecution at m — 2, then using the commutativity relations 
of the matrices we can write 


where (io, ii,..., im- 2 ) is a permutation of (0,1,..., m — 2). By (13.131) we can write 

= diag(/„,M^)Mm_i, (3.14) 

where t : {0,1,..., m — 2} — {1,..., m — 1} is a bijection such that for i = 0 : m — 3, 
T has a consecution (resp., inversion) at i if and only if a has a consecution (resp., 
inversion) at i. So by the induction hypothesis, Mr = Wm-s- Finally the product 
given in (13.1411 is 



—1 in 

■ 

0 

0 

0 


In 0 


0 lFm-3(:,l) Wm-3(:, 2 : m - 1) Wm-3(:,m) 


i (m —2)n 




Ir \ 


“_VF^_3(:,1) 0 VF™_3(:,2:m-l) VF™_3(:,m)J’ 

which is equal to Wm -2 constructed when cr has a consecution at m — 2. 

Case II. If a has an inversion at m — 2 the proof is similar, but 

= M™_idiag(/„, M^), 

since M^-i lies to the left of Mm_ 2 . □ 

4. Fielder pencils are linearizations for system polynomials. We now 

show that a Fiedler pencil Lo.(A) of the Rosenbrock system polynomial 5(A) is a 
linearization of 5(A) in the sense of system equivalence. An n-by-n matrix polynomial 
A'(A) is said to be unimodular if det(Ar(A)) is a nonzero constant for all A S C. Let 
Gi(A) and G 2 (A) be n-by-n rational matrix functions. Then Gi(A) is said to be 
unimodularly equivalent to G 2 (A) (denoted as G(A) G 2 (A)) if there are n-by-n 
unimodular matrix polynomials U{X) and F(A) such that U{X)Gi{X)V{X) = G 2 (A) 
for all A, see [l4] . 

Definition 4.1 (System polynomial). An {n + r) x (n + r) matrix polynomial 
X(X) is said to be a system polynomial of degree m if it is of the form 


\ XniX) 

X12 

X21 

X 22 (A) 


where Xu (A) is an n-by-n matrix polynomial of degree m, X 12 and X 21 are constant 
matrices and X 22 (A) is an r-by-r matrix pencil. If, in addition, the pencil X 22 (A) is 
regular then X(X) is said to be a Rosenbrock system polynomial and, in such a case, 
the rational matrix function 

r(A) := Xii(A) — X 12 (X22(A)) ^ X 21 

is said to be the transfer function associated with X{X). 

A system polynomial represents a linear time-invariant system in state-space form 
and is also referred to as a system matrix. Observe that if X{X) is an (n-l-r) x {n-\-r) 
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system polynomial then for any n-by-n unimodular matrix polynomials U{X) and 
V (A), the matrix polynomial 


Z(A) 


r u{x) 

0 

'VtW 

r v{x) 

0 

0 

Ir 


0 

Ir 


need not be a system polynomial. We now introduce the notion of a system equivalence 
which is different from “strict system equivalence” introduced by Rosenbrock [14) . 

Definition 4.2 (System equivalence). Let ft’(A) and iV(A) be (n + r) x (n + r) 
system polynomials. Then X{X) is said to be system equivalent to iV(A) (denoted as 
X{X) ~se y(.X)) if there are n-by-n unimodular matrix polynomials U{X) and R(A) 
such that for all A S C, we have 


U{X) 

0 

V( \\ 

r F(A) 

0 

0 

Ir 

A. 

0 

Ir 


y{x). 


Suppose that X{X) and 3^(A) are (n + r) x (n + r) system polynomials and that 
-^(A) ~se iV(A). Then it follows that X 22 {X) = A 22 (A) for A e C, where X 22 {X) and 
Y 22 {X) denote the (2,2) blocks of X{X) and 3^(A), respectively. Consequently, if X{X) 
is a Rosenbrock system polynomial then so is 3^(A). Let Tx{X) and Ty{X), respec¬ 
tively, denote the transfer functions associated with d:’(A) and 3^(A). Then it is easy 
to see that the transfer functions Tx{X) and Ty{X) are unimodularly equivalent. Con¬ 
sequently, system equivalence preserves (finite) transmission zeros (and their partial 
multiplicities) of a linear time-invariant system in state-space form. 

If T’(A) is an {n -\- r) x {n -\- r) system polynomial then for any positive integer p 
the (n -\- p r) x (n -\- p r) system polynomial given by 


Ip © X{X) := 


Ip 0 
0 X{X) 


is called an extended system polynomial m- This amounts to introducing p addi¬ 
tional zero output variables to the linear time-invariant system associated with X(X). 
We now introduce the notion of a system linearization for system polynomials which 
we refer to as Rosenbrock linearization. 

Definition 4.3 (Rosenbrock linearization). Let X{X) be an (n + r) x (n + r) sys¬ 
tem polynomial of degree m > 1 . Then an {nm-\-r) x (nm + r) system pencil L(A) (sys¬ 
tem polynomial of degree 1) is said to be a Rosenbrock linearization ofX{X) i/L(A) ~se 


I{m-i)n © X{X). In other words, L(A) is of the form L(A) = 
there are nm x nm unimodular matrix polynomials U{X) and R(A) such that 


L{X) 


Y 


X 


A-XE 


and 


\ U{X) 

0 

0 

Ir \ 


L(A) 


V{X) 

0 



0 

0 

Ir \ 


0 

X{X) \ 


where L{X) is an nm x nm matrix pencil, A — XE is an r-by-r matrix pencil and, X 
and Y are constant matrices. If, in addition, U{X) and V{X) are constant matrices 
then L(A) is said to be a strict Rosenbrock linearization of X{X). 

Observe that L(A) is a “trimmed structured linearization” of T’(A). Indeed, L(A) 
is a structured linearization in the sense that ]L(A) is itself a system matrix and hence 
is associated with a linear time-invariant system in SSF. Further, L(A) is a trimmed 
linearization of X(X) in the sense that L(A) is an (mn + r) x (mn + r) pencil whereas a 
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usual linearization of X{X) would result in a pencil of dimension m{n + r) x m(n + r). 
Furthermore, a byproduct of Rosenbrock linearization of X{X) is that it deflates some 
(all) of the infinite eigenvalues of X{X) when deg(ff(A)) > 2. Indeed, suppose that 
X{X) is regular and that deg(ff(A)) > 2. Then obviously T’(A) has eigenvalues at oo. 
On the other hand, L(A) has an infinite eigenvalue if and only if at least one of the 
pencils A — XE and L{X) has an infinite eigenvalue. Consequently, L(A) does not have 
an infinite eigenvalue when the leading coefficients of the matrix polynomial X{X) 
and E are nonsingular even though the system polynomial T’(A) has eigenvalues at 
oo. Thus L(A) deflates all infinite eigenvalues of X{X) when deg(fF(A)) > 2 and the 
leading coefficient of Al(A) is nonsingular. 

Let T’(A) be an (n + r) x (n + r) system polynomial of degree m. Suppose that 
L(A) is a Rosenbrock linearization of T’(A). Then T’(A) and L(A) are of the form 


X{X) = 


X{X) 

X12 

QTirl ir (W 

\ L{X) 

Li 2 

X21 

Y -XZ 


L21 

Y-XZ 


(4.1) 


where X{X) is an n-by-n matrix polynomial of degree to, L{X) is an mn x mn matrix 
pencil, Y — XZ is an r-by-r matrix pencil and ^ 12 ,^ 21 ).^i 2 jare matrices of 
appropriate dimensions. It is easy to see that L{X) is unimodularly equivalent to 
I(m-i)n®X{X), that is, L{X) is a linearization [71 [Tl] of the matrix polynomial X(A). 
Thus L{X) is a linearization of X{X) when L(A) is a Rosenbrock linearization of X{X). 

On the other hand, suppose that L{X) is a linearization of ^(A). Let t/(A) 
and R(A) be mn x mn unimodular matrix polynomials such that U{X)L{X)V{X) = 
I(m-i)n ® X{X) for all A S C. Then we have 


\ U{X) 

0 

0 

Ir \ 


L(A) 


■ ^(A) 

0 



0 

0 

Ir 


0 

X{X) 


(4.2) 


where L(A) 


\ m 

C/(A)-i(e^OAi2) 1 

_ (e^®A2i)H(A)-^ 

Y-XZ 


. Since a system poly¬ 


nomial is in general not invariant under system equivalence, it follows that L(A) 
is a Rosenbrock linearization of X{X) if and only if U{X)~^{em O -^ 12 ) and (e^ O 
-^ 2 i)I^(A)“^ are constant matrices independent of A. 


Corollary 4.4. Let X{X) and L(A) be as in B- If L(A) is a Rosenbrock 
linearization of X{X) then L{X) is a linearization of X{X). Conversely, if L{X) is a 
linearization of X{X) then L(A) is a Rosenbrock linearization of X{X) if and only if 
U{X)~^{em ^ X 12 ) and {ej^ ^ X 2 i)V{X)~^ are constant matrices (independent of X) 
and are equal to L 12 and L 21 , respectively. 

So, the moot question is: Are U{X)~^{cm® X 12 ) and (e^ O A 2 i)R(A)“^ constant 
matrix polynomials when L{X) is a linearization of A(A)? It turns out that it is 
indeed the case when L{X) is a Fiedler pencil of A(A) and that the Fiedler pencils 
in Theorem 1,4.61 are Rosenbrock linearizations of A)A). Computation of the last block 
column and row of U{X)~^ and I4(A)“^, respectively, is a laborious task and involves 
long calculations. We therefore adapt the method of proof presented in |5] for matrix 
polynomials to the case of system polynomials and show that Fiedler pencils of system 
polynomials are indeed Rosenbrock linearizations. 

Definition 4.5 (Horner shift, [6]). Let P{X) = Aq + AAi -I- • • • -I- X'^Am be a 
matrix polynomial of degree to. For fc = 0,... to, the degree k Horner shift of P{X) is 
the matrix polynomial Pk{X) := Am-k + XAm-k+i -f ■ • • -f X^Am- These Horner shifts 
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satisfy the following: 


Po{^) — Ajn, 

Pfe+i(A) = APfe(A) + Am-k-i, for 0 < k < m - 1, (4.3) 

P^(A) = P(A). 


Let H := {Hij) be a block m x n matrix with p x q blocks Hij. The block 
transpose of H, denoted by , is the block n x m matrix with p x q blocks defined 
by := Hji, see [6]. Consider the auxiliary matrices associated with a matrix 

polynomial given below. 

m 

Definition 4.6 (i, Definition 4.2). Let P{\) = 'Yf, AMi be an n x n matrix 

2—0 

polynomial, and let T’i(A) be the degree i Homer shift of P(X). For 1 < i < m — 1, 
define the following nm x nm matrix polynomials: 


Qi(A) 


I(i—l)n 

In XI 
On In 

n 

, RiiX) ;= 

On 

Ir. 

In 

P{X) 

- 

^{m — i—l)n 



-^(m — 2 —l)n 

Ti{X) := 

0(2—l)n 

0^ 

A7n 

APi-i(A) 

X^Pi-i{X) 


1 




0(m — 2 —l)n 



A (A) 


0(i —l)n 

L’i-l(A) On 

On In 

I{m — i — l)n 


and Dm[X) := diag [0(m-i)n, Ai-i(A)] . For simplicity, we often write Qi, Ri,Ti, Di 
in place of Qi{X),Ri{X),Ti{X),Di{X). Note that Di{X) = Mm, and Qi{X), Ri{X) are 
unimodular for all i = 1,. .. ,m — 1. Also note that Rf{X) = Ri(X). 

The auxiliary matrices satisfy the following relations. 

Lemma 4.7 m, Lemma 4.3). Let Qi, Ri,Ti, Di be as in Definition and 
Mi’s be Fiedler matrices associated with P{X). Then the following relations hold for 
i = 1 ,.. . ,m — 1 . 

(^) Qf i^X)i)Ri = XDi^i + Ti, and Qf {Mm-(i+i)Mm-i)Ri = Mm-(i+i) + Ti. 

(b) Rf {XDi)Qi = XDi^i + T®, and Rf{Mm-iMm-{i+i))Qi = Mm-(i+i) + T^. 

(c) TiMj = MjTi = T, and T^Mj = MjT^ = for all j <m-i-2. 

We extend these results to the case of system polynomials. For this purpose, we 
introduce “block transpose” of a system matrix as follows. 

Definition 4.8 (Block transpose). Let A be an {mn + r) x {mn + r) system 
matrix given by 


xl:= 


A 

Ci® X 

ej (g) Y 

Z 


where A := [Aij\ is an m x m bloek matrix with Aij € C"^", X € € C” 


Z G and et is the k-th column of I„ 

is defined by 


The block transpose of A, denoted by A^ 


\ 

ej 0 X 

ef ® Y 

Z 
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where is the block transpose of A. 

Recall that C'p(A) is the first companion form of the matrix polynomial -P(A). It 
is well known that the second companion form (72 (A) of P{X) is the block transpose 
of Cp{X), that is, (72(A) = (7p(A)®, see [6]. It turns out that the same is true for 
the second companion form C 2 (A) of the system polynomial S{X). Indeed, in view of 
Theorem 13.61 Proposition 13.71 and (I3.12L we have C 2 (A) = C(A)®. 

We now define auxiliary system polynomials as follows. 

Definition 4.9 (Auxiliary system polynomials). Let Qi{X),Ri{X),Ti{X), and 
Di{X) be as in Definition \4-6\ Fori = 1,... ,m—1, define (nm + r) x (nm + r) system 
polynomials: 


and 


Q^{X) := 
TfiX) := 
VmiX) := 


■ Q^{X) 

0 

0 


Ir \ 

r ux) 

0 

0 


Or 

' Dm 


0 

0 


-E 


^l^{X) 


V^{X) 


■ RfiX) 

0 


0 

r A (A) 

Ir \ 
0 

1 

0 

-E 



Note that T’i(A) = 


r A (A) 

0 


' Mm 

0 

0 

-E 


0 

-E 


and that Qi(A) and 


TZi{X) are unimodular matrix polynomials for i = I,...,to — I. Also, note that 
TZf{X) = TZi{X) for i = 1,... ,m — 1. 

The auxiliary system polynomials satisfy the following relations. 

Lemma 4.10. Let Qi,TZi,Ti,'Di be the system polynomials given in Definition 
\4-9\ and ’s be Fiedler matrices associated with S{X). Then the following system 
equivalence relations hold for i = 1 ,... ,m — 1 . 

(a) QfiXVfiTZi = XDi+i +Ti, and +71. 

(b) 7^“(AA)2^ = AP.+i + T;®, and 7e®(M™_,M,„_(,+i))Q, = + if. 

(c) TMj = MjTi = % and = M^T)® = 7)® for all j <m-i- 2 . 

Proof (a) We have 


Q®(AI?,)7^, 


Qf 

0 

\ aa 

0 

Ri 

0 


\ QfXD,R, 

0 

0 

Ir 

0 

-XE 

0 

Ir 


0 

-XE 


XDi^i + Ti 


0 


-XE 


— XDi-\-i + 71 


and 


Qf(Mm_(i+i)Mm,-i)7?.i — 


\ Qf 

0 


0 

' Mm-^ 

0 

Ri 

0 

0 

Ir 

0 

Ir 

0 

Ir 

0 

Ir 


Qf Mm-(i+l)Mm-iRi 

0 


Mm-(i+l) + Ti 

0 

0 

Ir 


0 

Ir 


Mm-(z+i) + Ti. For 
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i = m — 1, we have Q 5 ^_ 2 (MoMi) 7 ?.m_i 



Similarly 7^®(Mm_iM„_(J+l))QJ = M„_(i+i) + T"®. 
(c) Finally, we have 



Similarly 7'®Mj- = = 7^®. For i = m - 2 we have 



completes the proof. 0 
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As in the case of matrix polynomial [B] , we now define a family of system pencils 
which will form the intermediate steps in the system equivalence transformation of 
a Fiedler pencil Lcr(A) into I(m-i)n ©>S(A). Note that a Fiedler pencil of S{X) is a 
system pencil. 

Definition 4.11. Letha-iX) = AM^— be the Fiedler pencil ofS{X) associated 
with a hijection a. For j = 1, 2,..., m, define 

M« := n 

(z)<m—J 


where the factors Mcr-i(i) are in the same relative order as they are in Mo-. Note that 
= no-i(i)<m-i thatM.’^T'^ = Mq. Also for j = l,2,...,m 

define the (nm + r) x (nm + r) system pencils Lct^(A) := X'Dj{X) — Observe that 

Lct^^(A) = XVi — = AM„i — Mo- = Lcr and that 


lM(A) = AP, 


Dm 

0 


^(7n — l)n 

0 

0 

-E 


0 

5(A) 


The next result shows that (A) ~se (A) for f = 1, 2,..., m — 1. 

Lemma 4.12. We have Lct^(A) ~se Lct^^^(A) for i = 1,2,. ..,m — 1. More 
precisely, if Qi and TZi are the system polynomials given in Definition \4-9[ then 


L(?+i)(A) 


QfL«(A)7^„ 


if a has a consecution at m — i — 1 
if a has an inversion at m — i — 1. 


Proof. The proof is exactly the same as that of Lemma 4.5 in i- □ 

It is now immediate that a Fiedler pencil is a Rosenbrock linearization of >S(A). 
Theorem 4.13 (Rosenbrock linearization). Let 5(A) be an (n + r) x (n + r) 
system polynomial (regular or singular) of degree m. Then a Fiedler pencil Lct(A) of 
the system polynomial 5(A) is a Rosenbrock linearization of S{X). 

Proof. By Lemma 14.121 we have m — 1 system equivalences 


L.(A)=L«(A) Lf (A) 


lW(A) 



0 

0 

5(A) J 


(4.4) 


where Lct^(A) is as in Lemma [4. 121 This shows that Lct(A) ~se I{m-l)n ® 5(A). □ 

The m — 1 system equivalences in (14.41) provide the desired system equivalence 
that transforms L£r(A) to the extended system polynomial I(m-i)n ©5(A) and can be 
constructed as follows. 

Corollary 4.14. Let Lo.(A) be the Fiedler pencil o/ 5(A) associated with a 
bijection a, and Qi, TZi for i = 1, 2,... m — 1, be as in Definition\4.S\ Then 


U{X)U{X)V{X) = 


^(m — l)n 

0 

0 

5(A) J 


(4.5) 


where U{X) and V(A) are {nm + r) x (nm + r) unimodular system polynomials given 
by 


a 


'm—(i+l) ’ 
^-( 2 + 1 )’ 


U{X) :=U[)Ui-■ ■Um-d.Um-2, withUi = 
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if a has a consecution at i, 
if a has an inversion at i. 






















V(A) := V,n-2Vm-3 ' ' ' VlVo, with H = 


if a has a consecution at i 
if a has an inversion at i. 


(i+l)) 

The indexing oiUi and Vi factors inZ^(A) and V’(A), respectively, in Corollarv l4.14l 
has been chosen for simplification of notation and has no other special significance. 

Example 4.15. For m = 4 and Lct(A) = AM 4 —M 2 M 0 M 1 M 3 , by Corollary\4.14\ 
we have (Q 57 ^ 2 Q?)L,(A)( 7 ^lQ 27 ^ 3 ) = -Isn ©5(A). ■ 

Observe that the invariant zeros of an LTI system in SSF such as E in dm can be 
computed in two steps. First, construct a Fiedler pencil Lcr(A) using Algorithm[TJ Sec¬ 
ond, solve the GEP Lcr(A)it = 0 by standard algorithm such as the QZ-algorithm. On 
the other hand, the decoupling zeros of the LTI system can be obtained by computing 
eigenvalues of the pencils in (12.2p using standard algorithm such as GUPTRI [H [^ . As 
for the transmission zeros of the LTI system, it is well known m that they coincide 
with the invariant zeros when the LTI system is controllable and observable, see next 
section. Thus we have a direct method for solving Problem-I. 

5. Linearizations for rational matrix functions. We now introduce “lin¬ 
earization” of a rational matrix function for solving a rational eigenvalue problem. 
Intuitively, a linearization of a rational matrix function G(A) should be a matrix pen¬ 
cil L(A) of appropriate dimension (ideally the smallest dimension) such that G(A) and 
L{X) have identical “zero structure” in C, that is, they have the same spectrum and 
that the partial multiplicities of the zeros (eigenvalues and eigenpoles) of G(A) are the 
same as those of L{X). In other words, Sp(G) = Sp(L) and Ind 0 (A,L) = Ind 0 (A, G) 
for each A € Sp(L). Additionally, a linearization L{X) should allow an easy recovery of 
eigenvectors of G(A) from those of L{X). To achieve this goal, we proceed as follows. 

Let G(A) € C(A)’”^" be a rational matrix function. Then G(A) is said to be 
proper if G(A) —> Gq G as A —>• 00 . On the other hand, G(A) is said to be 

strictly proper if G(A) —>■ 0 as A ^ 00 , see [T31[T1]. A rational matrix function G(A) 
can be uniquely decomposed as 


G(A) = P(A) + Q(A), 


(5.1) 


where P{X) is a matrix polynomial and Q(A) is a strictly proper rational matrix 
function. We denote the degree of P(A), the highest power of A in P(A), by deg(P). 
Recall that SM(G(A)) denotes the Smith-McMillan form of a rational matrix function 
G(A) and SF(P(A)) denotes the Smith form of a matrix polynomial P{X). Also recall 
that ifaW is the pole polynomial of G(A). The degree of the pole polynomial V'g(A) 
is called the least order of the rational matrix function G(A), see [T5] . 

Definition 5.1 (Linearization for rational matrix function). LetG{X) G C(A)"^"' 
he a rational matrix function such that G(A) = P(A) + Q{X), where P(A) is a matrix 
polynomial and Q{X) is a strictly proper rational matrix function. Suppose that 


SM(G(A)) = diag 


<(>i(A) 02 (A) 0 fc(A) 
'01 (A)’ 02 (A)’'"’ 0 fe(A) 


,0, 


where k = nrk(G). Set m := max(deg(P), 1) and r := deg(0G)j where 0g(A) is the 
pole polynomial of G(A). Then an (mn -I- r) x fmn -\- r) matrix pencil L(A) = XX + y 
is said to be a linearization of G{X) if the Smith form o/L(A) is given by 


SF(L(A)) = diag 0i(A), • • • ,0fc(A),O,... ,0) . 
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We show that a rational matrix function admits a linearization and that such a 
linearization can be constructed from a minimal state-space realization of the rational 
matrix function. A state-space realization of a rational matrix function G(A) is an 
LTI system in SSF for which G(A) is the transfer function [HI IT^ [21 ITU]. 

Definition 5.2 (State-space realization). A tuple S := (AG — A, B, C, G(A)) G 
C[A]''^’’xC'’^" xC[A]”^^” with E being nonsingular is said to be a realization 

of a rational matrix function G(A) S C(A)"*^" if G(A) = P(A) -I- C{XE — A)~^B. In 
such a case, G(A) is said to be the transfer function of the realization S. The matrix 


polynomial S{X) := 


■ PW 

G 

B 

-< 

1 


is said to be the Rosenbrock system 


polynomial associated with the realization S. The LTI system in SSF associated 
with the realization S, which we also denote by E, is given by 


Ex{t) = Ax{t) -\- Bu{f) 
y{t) = Cxit) + P{f^)u{t) 

for which G(A) is the transfer function andS{X) is the Rosenbrock system polynomial. 

We mention that the nonsingular matrix E is can be chosen to be the identity 
matrix A. Since G(A) = f’(A) -|- Q{X) and Q{X) is a strictly proper, there exist r x r 
matrices A and E with E being nonsingular, B G and G G such that 

Q(A) = C{XE - A)-^B, see [11]. Thus G(A) = P(A) -b C{XE - A)-^B. Obviously, 
a state-space realization of G(A) is not unique. We refer to [U [TOl [H] [19] for theory 
and algorithms on realizations of rational matrix functions. 

A realization (AG — A, B, G, G(A)) of G(A) is said to be minimal if the size of the 
pencil AG — A is the smallest [14] [19] . The minimality of a realization is characterized 
by controllability and observability of the associated LTI system in SSF. 

Theorem 5.3. Let T. := (AG - A, B, G, G(A)) € C[A]’''<'’ x C’'"'" x 

C™x’’xC[A]’”X" a realization ofG{X). Then the following conditions are equivalent. 

(a) The realization E is minimal. 

(b) The LTI system associated with the realization E is controllable and observ¬ 
able. Equivalently, 

rank ([AG — A G]) = r and rank 

(c) The size of XE — A is equal to the least order o/G(A), that is, r = deg{ipG), 
where ifciX) is the pole polynomial ofG{X). 

A realization of G(A) which is not minimal can always be reduced to a minimal 
realization m Ho]. Observe that if (AG — A, B, G, G(A)) is a minimal realization 
of G(A) then the associated LTI system E does not have input decoupling zeros as 
well as output decoupling zeros. In such a case, A S C is a pole of G(A) if and only 
if A is an eigenvalue of the pencil AG — A. Further, the zeros of G(A) coincide with 
the eigenvalues of the Rosenbrock system polynomial 5(A) associated with the LTI 
system E. In other words, the transmission zeros of the LTI system E coincide with 
the invariant zeros of the system. This can be seen as follows. 

Theorem 5.4. Let P(A) be an (r -|- m) x {r -\-n) matrix polynomial given by 


XE-A 
C 


= r for all X G C. 


P(A) := 


r -T(A) 

U{X) 1 

[ V{X) 

W{X) _ 
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--nx) 
[ . 

for all A € C. Consider the m x n rational matrix function G(A) given by 


where T(A) € C[A]'’^’’ is regular and rank([—T(A) C/(A)]) = rank 


G(A) := W{X) + V{\){T[X))-^U{\). 


Suppose that the Smith-McMillan form of G{X) is given by 

SM(G(A)) = diag ((/.i(A)M(A),..., <^,(A)/V^,(A), . 

Then the Smith forms ofT{X) and P(A) are given by 

SF(r(A)) = Aia.g{Ir-q,lfq{X),lfq-l{X) . . . ,'lfl{X)) , 
SF(P(A)) = diag (4, (/)1 (A), 02 (A) . . . ,(j)q{X),Om-q,n-q) ■ 


For a minimal realization of a rational matrix function we have the following. 
Theorem 5.5. Let {XE - A, B, G, P(A)) G C[A]’'^’' x C''"'" x x C[A]™''" 
be a minimal realization of G{X) G C(A)'"^"' and 5(A) be the associated Rosenbrock 
system polynomial. Consider the pencil L{X) := XE — A and suppose that 


SM(G(A)) = diag 


01 (^) 02(A) 0fc(A) 

0i(A)’02(A)’-'-’0fe(A)’ , 


(a) Then the Smith form SF(L(A)) of L{X) is given by 


SF(L(A)) = diag {Ir-k, 


and the Smith form SF(5(A)) o/5(A) is given by 


SF(5(A)) = diag(4,01,02,...,0fe,0,...,0). 

(b) We have Poles(G) = Sp(L) and Ind^(A, G) = Ind,^(A, L) for X G Sp(L). 

(c) We have Sp(G) = Sp(5) and Ind,^(A, G) = Ind,^(A,5) for X G Sp(5). 

(d) fTe have Eip(G) = Sp(5) n Sp(L) and Eig(G) = Sp(5) n (C \ Sp(L)). 

Proof. Since the realization is minimal, the results in part (a) follow from Theo¬ 
rem 15.31 and Theorem 15.41 On the other hand, the results in parts (b), (c) and (d) 
follow from part (a) and Theorem 12.31 □ 

We now show that Fiedler pencils of the Rosenbrock system polynomial associated 
with a minimal realization of G(A) are in fact linearizations of G(A). 

Theorem 5.6 (Linearization). Let {XE - A, B, G, P(A)) G C[A]’''''’ x x 
C"^’"xC[A]"'^" be a minimal realization ofG{X) G C(A)"'^" andS{X) be the associated 
Rosenbrock system polynomial. Suppose thatm is the degree of P{X). TetLcr(A) be the 
Fiedler pencil ofS{X) associated with a bijection ct : {0,1,..., m — 1} {1, 2,..., m}. 

Thenh^{X) is a linearization of the rational matrix function G(A). 

Proof. Suppose that the Smith form of G(A) is given by 


SM(G(A)) = diag 


01 (A) 02(A) 0fc(A) 
01 (A)’ 02(A)’'"’ 0fe(A) 


,0, 


By Theorem 15.51 the Smith form SF(5(A)) of 5(A) is given by 


SF(5(A)) = diag (1^,01,02,..., 0fc, 0,..., 0). 
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Since L£r(A) is a Rosenbrock linearization of 5(A), there are nm x nm unimodular 
matrix polynomials 17(A) and R(A) such that 


r u{x) 

0 

0 

Ir \ 


U(A) 


17(A) 

0 



0 

0 

Ir 


0 

5(A) 


As 5(A) is unimodularly equivalent to SF(5(A)), it follows that Lcr(A) is unimod- 
ularly equivalent to diag (/(m-i)ra+r) ^^* 1 ) 02) ■ • ■ > 4>k, 0,..., O) . Hence the Smith form 
SF(L<^(A)) is given by 

SF(L^(A)) = diag 0i, 02, ■ ■ ■, 0fc, 0,..., O) 


proving that Lcr(A) is a linearization of G'(A). □ 

Thus we conclude that the zeros of a rational matrix function G(A) G C(A)"'^" 
in general and the eigenvalues of G(A) in particular can be computed by solving 
the generalized eigenvalue problem for a Fiedler pencil of the Rosenbrock system 
polynomial associated with a minimal realization of G(A). Hence a direct method for 
solving a rational eigenvalue problem G{X)u = 0 can be summed up as follows. 


Algorithm 2 A direct method for solving an REP G{X)u = 0. 

1. Compute a minimal realization E := {XE — A, B, G, P{X)) of G(A). 

2. Compute a Fiedler pencil Lcr(A) of the Rosenbrock system polynomial asso¬ 
ciated with S using Algorithm [T] 

3. Solve the generalized eigenvalue problem Lcr(A)u = 0. 


Note that the eigenvalues of Lct(A) are precisely the eigenvalues and eigenpoles 
of G(A), that is, Sp(G) = Sp(Lct). Also note that the eigenpoles of G(A) are those 
eigenvalues of Lct(A) which are also eigenvalues of the pencil L{X) := A — XE, that is, 
Eip(G) = Sp(Lcr) n Sp(L). We mention that the pencil Lo-(A) allows easy recovery of 
eigenvectors of G(A) from those of Lo-(A) which is considered in a forthcoming paper. 

6. Conclusions and future work. Computation of zeros of a linear time- 
invariant (LTI) system is a major task in Linear Systems Theory. Eor computing 
zeros of an LTI system E in SSE, we have introduced Rosenbrock linearizations, 
which are trimmed structured linearizations, of the Rosenbrock system polynomial 
5(A) associated with the LTI system E. We have introduced Fiedler-like pencils for 
the Rosenbrock system polynomial 5(A) and have shown that the Fiedler-like pencils 
are in fact Rosenbrock linearizations of 5(A). For solving rational eigenvalue prob¬ 
lems, we have introduced linearizations of a rational matrix function by utilizing a 
minimal state-space realization of the rational matrix function. We have introduced 
Fiedler-like pencils for a rational matrix function by embedding the rational eigen- 
problem into an LTI system in SSF which is controllable as well as observable and 
have shown that the Fiedler-like pencils are in fact linearizations of the rational matrix 
function. Other important issues such as the eigenvector recovery of a rational ma¬ 
trix function from that of its linearizations, structured preserving Fiedler-like pencils 
(such as generalized Fiedler pencils and generalized Fiedler pencils with repetition) 
for Rosenbrock system polynomials and rational matrix functions will be considered 
in forthcoming papers. 
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